Lattice QCD calculation of K^K scattering length 
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We deliver ab initio calculation of s-wave K'^K~ scattering length {a^^^ ) by Liischer's formula. 
In the "Asqtad" improved staggered dynamical fermion formulation, we measure K^K~ four-point 
correlation function by moving wall sources without gauge fixing, and find ^ = 0.456±0.272 fm, 
which is in reasonable agreement with tree-level prediction and comparable with experimental result. 
An essential ingredient in our calculation is to explicitly include the disconnected diagram. 
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PACS numbers: 12.38.Gc 

K^K^ interactions are still very poorly understood. 
Our basic motivation for investigating K'^K~ scatter- 
ing via lattice QCD is preliminarily discussed in Ref. [ij, 
as an attempt to unveil the mysterious nature of the 
scalar resonances ao(980)//o(980), as well as the (/)(1020), 
whose masses are pretty close to the K^K~ threshold. 
Some experiments of K'^K" scattering are conducted at 
the cooler synchrotron COSY facilities 0-01' where the 
prompt K^K~ spectrum is expressed by a four-body 
phase-space distribution altered by various final state in- 
teractions, which indicates a strong evidence of K~^K^ 
background. Silarski et al. first extracted K^K~ scat- 
tering length: a^+if- = [(0.2t|];^) + ^(0.41°:!^)] fm 0, 
via the low energy invariant mass distributions 

and the generalized Dalitz plot analysis 0, 01- Moreover, 
Ablikim et al. consider p meson decay to K'^K~ if its 
mass is larger than /^+A'~ threshold [8|. A clear signal in 
the /2(1525) region as well as a nonresonant component 
was observed in the KK mass spectrum A partial 
wave analysis of the centrally produced K~^K~ system 
has been carried out to observe the peaks in the s-wave 
corresponding to the /o(1500) and /o(1710), even the d- 
wave indicating the evidence for the /2(1270)/a2(1320), 
/;^(1525) and /2(2150) The possibility to measure 

Xco in the K^K~ channel are recently reported in 
Furthermore, KK interactions are very fundamental to 
many other interesting physics phenomena, e.g., a study 
of neutron star properties due to possible Bosc-Einstcin 
condensation |12| since kaons are bosons, and there was 
a pioneering lattice study on kaon condensation [l3l |. 

Although the study of K^K^ scattering length from 
QCD is a non-perturbative problem in essence, some the- 
oretical efforts are stiU being made to study KK/K~^K~ 
scattering [14 To understand the intrinsic nature of 
scalar mesons in terms of quark and gluon constituents 
the K~^K~ couplings of the isoscalar scalar mesons cr and 
/o(980) are investigated [l7| . 

Lattice QCD is at present the unique non-perturbative 
way which allows us to evaluate K~^K~ scattering length 
from first principles. Nevertheless, so far, no lattice QCD 
study on it has been reported because it is very difficult 
to compute vacuum diagram |H. Motivated by NPLQCD 



benchmark work on K^K^ scattering [iSf and our ex- 
ploratory work on KK scattering in the 1=1 channel [l[ , 
we further probe K~^K~ scattering from lattice QCD. 
The computations are launched on a 20"^ x 48 MILC full 
QCD gauge configuration with the 2 + 1 flavors of the 
Asqtad improved staggered sea quarks at TOtt w 240 Mev 
and lattice spacing a w 0.15 fm |19l.[20j. 

Let us learn scattering of two Nambu-Goldstone kaons 
with zero momentum in the staggered fermion formal- 
ism pH - fisf . Using operators Ok+{xi) and Ok-{x2) 
for if+ and K~ at points xi, X2, respectively, we write 
K'^K~ four-point function as 

^ (Ok+(x4)Ok- (-^3)0],+ (2^2)0],- (2:1)), 



Xi,X2 X3,X4 



with the kaon interpolating field operators defined by 
^K+ix) — s{x)j^u{x) and Oj^-(t) = u{x)j5s(x), where 
xi = (xi,ti), and likewise for X2, X3 and X4. In the 
isospin limit, the diagrams contributing to the four-point 
function Cj^+j^- corresponding to the terms (l)-(4) in 
Appendix A of Ref. [l[ are shown in Fig. 1 of Ref. M . To 
prevent the Fierz rearrangement of quark lines [2^ [23j , 
we opt ti = 0, ^2 = 1, ^3 = t, and + In Fig. [1] we 

illustrated the quark line diagrams, which are labeled as 
direct (£>), rectangular and vacuum (V) dia- 

grams, respectively [l|. The direct diagram can be easily 
calculated. However, other three diagrams require extra 
quark propagator connecting time slices ^3 and [l| . 
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(a) D (b) i?„ (c) R, (d) V 

FIG. 1: (color online). Quark-line diagrams contributing to 
K'^ K~ four-point functions. Short bars show wall sources. 
Open circles are sinks for local K'^ / K~ operators. Blue and 
red lines are u/d and strange quark lines, respectively. 
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We compute these diagrams by the moving wall sources 
technique introduced by Kuramashi et al [22|, |23| . The 
concrete lattice calculation of these diagrams are de- 
scribed in Ref. [l[ . The rectangular and vacuum diagrams 
create gauge- variant noise [22, 23|, and we usually reduce 
it by conducting gauge field average without gauge fix- 
ing jll. l25l - l28j . In terms of four diagrams in Fig.[Tl K^K^ 
correlation function can be expressed as 

Ck+ k- {t)=D- NfRu - NfR, + V, 

where the staggered-flavor factor Nf is inserted to correct 
for the extra factors Nf in valence fermion loops [2l| . The 
four-fold degeneracy of staggered sea quarks is removed 
by taking the fourth root of the staggered fermion deter- 
minant [21I, [2^. The s-wave K^K~ scattering length in 
the continuum is defined in terms of phase shift by [l| 

tan(5o(fc) . . 

ao ~ lim ; , where the scattermg is purely clas- 

tic below inelastic thresholds. Wc should remember that 
the truncation of the effective range r serves as a source 
of the systematic uncertainty which appears as 0{1/L^) 
and hence we ignore it [ij. The scattering momentum k 
is related to total energy of K~^K^ system by 



2. ml + k^ 



Er+K '"-K 

and (5o(fc) is evaluated by Liischer's formula [3C 



tan (5o(fc) 
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here the zcta function Zqo is formally denoted by 
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and can be calculated by the approach described in [3l| . 
The energy Ej^+j^- can be obtained from K^K~ four- 
point function with large t, whose behavior is given as 32 1 



Cx+K-it) = Zk+k- cosh 
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The oscillating term corresponding to the contributions 
from the intermediate states with opposite parity is a 
characteristic of the Kogut-Susskind formulation [3^ ■ 
In practice, we also calculate the ratios [l[ 



where Ck is kaon two-point function. We can verify that 
when t <^ T/2, even we impose periodic boundary con- 
dition in the temporal direction, we still can roughly ex- 
tract the energy shift from these ratios. 



We should bear in mind that the contributions of non- 
Nambu-Goldstone kaons in the intermediate states is ex- 
ponentially suppressed for large time region due to their 
heavier masses compared to these of Goldstone kaon [2l| . 
Therefore we neglect this systematic error. 

In this work we measured K~ four-point function 
on a MILC lattice ensemble of 601 20'^ x 48 gauge config- 
urations in the presence of the 2-1-1 dynamical flavors of 
the Asqtad-improvcd staggered fermions with bare quark 
masses ara^j^j anis = 0.00484/0.0484 and bare gauge 
coupling 10 /g^ = 6.566 [l9l.|20l. The inverse lattice spac- 
ing a" 



1.373;?^ GeV [19|,12 



The strange quark mass 
is fixed to its physical value [19|, [2^ . Periodic boundary 
condition is applied to spatial and temporal directions. 

We employ standard conjugate gradient method to 
compute the required matrix element of inverse fermion 
matrix. We compute the correlators on each of the 48 
time slices, and explicitly combine these results. After 
averaging the correlator over all 48 possible values, the 
statistics are greatly improved. For each time slice, six 
fermion matrix inversions arc required corresponding to 
the possible 3 color choices for each kaon source. Totally 
we conduct 288 inversions on a single configuration. 

In Fig. m the individual ratios, {X = D,Ru,Rs 
and V) are shown as the functions of t. The value of 
the direct amplitude R^ is pretty close to unity. After 
an initial increase up to i ~ 4, the rectangular R^" and 
R.^'" amplitudes demonstrates a roughly linear decrease 
up to i '-^ 18, and the oscillating behavior in large time 
region is clearly observed [s^ . The features are what we 
expected P, [U [s^l • We also notice that R^" oscillate 
strongly than i?^" . The physical meaning of the intense 
oscillation of R^" and its mass dependence is not clear 
to us, which are needed to be further investigated. 

The vacuum amplitude is pretty small up to t ~ 8 — 13, 
and loss of signals after that. Since its correlation func- 
tions reflect the quantum fluctuations of QCD vacuum, 
the error for disconnected amplitude should be roughly 
independent of t jsl], thus the error of R^ (t) grows ex- 
ponentially as e^"*^* [l| for large t. While the analytical 
arguments in 33| suggests the error of ratio for rectangu- 
lar Ru diagram grow exponentially as g(2™-ff-'Tis5)t^ where 
russ is a flctitious meson with two valence quarks mass 
about nis (likewise for Rg diagram as gi'^'^K-m^)t^ 1 j,^^ 
magnitudes of errors are in quantitative agreement with 
these predictions as exhibited in Fig. |3l Fitting these er- 
rors SR-^{t) by a single exponential SR^{t) ~ e^^* for 
8 < t < 18, we get afix = 0.246, 0.581 and 0.755 for 
X — Ru, Rs and V, which can be compared with our 
determinations of am^ = 0.1750(2), arriK = 0.3774(2) 
and aniss = 0.5012(2) on the same set of gauge configu- 



^ For the MILC gauge configuration used in we can easily 
check that m-n- 2m, — mss and rusa ~ 2m — m,r. 
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FIG. 2: (color online). Individual ratios it) {X = D, R^, 
Rs and V) for K~ four-point function as functions of t. 
(a) Direct diagram shifted by 0.5 (red diamonds) and vacuum 
diagram (blue octagons); (b) Rectangular diagram Ru (blue 
octagons); (c) Rectangular diagram Rs (red octagons). 



rations. In principle, we can reasonably assume that the 
vacuum amplitude is still small for large t, and ignore 



it in the following analysis (22|. Nervelessness, we will 



include it such that the calculation has integrity. 




FIG. 3: (color online). Error of ratios R {t) {X = Ru, Rs 
and V) as functions of t. Solid lines are single exponential fits 
over 8 < t < 16. 



compared with lattice data. As we expected, the sys- 
tematically oscillating feature in the large time region is 
clearly observed [s^. The scattering momentum k'^ in 
lattice unit calculated by Eq. ([T]) and then s-wave scat- 
tering length a^*^ obtained through Eq. (I2|), together 
with fitted values of aE, fitting range and fitting quality 
(x'^/dof) are listed in Table HI The errors of scattering 
momentum k and scattering length are calculated from 
the statistic errors of the energy aE and kaon mass amx- 



The errors quoted for a, 



K+R- 



are statistical only. 
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FIG. 4: (color online). The K'^ K~ correlation function. The 
fitting range is indicated by points and fitted lines in red 
(darker points and lines), (a) The result excluding vacuum di- 
agram, (b) The result including vacuum diagram. Occasional 
points with negative central values are not plotted. 

TABLE I: Summary of the s-wave K'^ K~ scattering length. 
arriK = 0.3777(2) and a/^ = 0.1177(4) are used to estimate 
the tree- level prediction value in Eq. ((5]). The energy E and 
scattering momentum k axe listed in lattice units. Column 
two show the results left vacuum diagram out. Column three 
display the results with the inclusion of vacuum diagram. Fit- 
ting range and fit quality (x^/dof) are listed as well. 

Lat. Results(Vout) Lat. Results(Vin) Tree Level 



Fit A 



Fit B 



K+K~ 

k" 

Range 
xVdof 



I. 122(170) 
0.74698(85) 
-0.00608(62) 
7-16 

II. 5/6 



0.873(520) 

0.74839(306) 

-0.00519(216) 

6-16 

9.0/7 



0.8188(61) 



We use Eq. to extract the energy aE, then insert it 
into Eq. ([2]) to obtain the scattering length. In practice, 
we select the energy aE by looking for a combination of 
a "plateau" in the energy as a function of the minimum 
distance Dmin and a good confidence level (x^/dof) for 
the fit d. The K+K' correlation function is plotted 
in Fig. m where we can watch fitted functional form as 



In Table U we also list the s-wave K^K scattering 



length predicted by the tree level xPT formula 



scatte ring 
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(5) 



where we substituted the corresponding kaon masse ttik 
and the pion decay constant /^r determined on the same 
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set of gauge configurations. We noted that the simulation 
results are consistent with Eq. ([5]) within 1 2 standard 
deviations, which we find pretty satisfactory considering 
the systematic uncertainties discussed above 

It is definitely surprising that the extracted scattering 
lengths are in reasonable agreement, within uncertain- 
ties, with the tree-level prediction of SU(3) xPT, consid- 
ering one open question that tree-level xPT calculation 
for K'^K^ scattering length are certainly not accurate 
enough since one has the /o and oq resonances close by 
which make the nonperturbative calculation imprecise, 
maybe one-loop calculation and potential contribu- 
tions from higher order effects are probably important. 
This is somewhat similar to the cleaner case of K'^K'^ 
scattering length [l^ where there are no resonances near 
by. However, we should keep firmly in mind that there 
still exist a large statistical error for Fit A (without in- 
cluding vacumm term) and big error in Fit B (including 
vacuum diagram) as we predicted in [l| . 

In summary, in this work, we conducted the realistic ab 
initio lattice QCD calculation of s-wave K^K~ scatter- 
ing length, where the rectangular and vacuum diagrams 
play a crucial role. We employed the technique of the 
moving wall sources without gauge fixing [22, 23 1 to cal- 
culate all the four diagrams, and observed a clear signal 
of the attraction. We, first from lattice QCD, obtained 
its s-wave scattering length rriKaQ ^ = 0.873(520), or 



with the tree level xPT predictions [H, Hi, 15 1. 

Moreover, M. Silarski et al. carried out a special fit for 
us and achieved the experimental value of the s-wave 
K+K 



341, 



scattering length: Uq ' ^ = 3.5l2'o 
which was obtained by setting the imaginary part of 
K~^K~ scattering length to be zero when fitting their 
experimental data 0, But due to the low statistics, 
the uncertainties are still rather large. It is quite encour- 
aging that with the current computing resources, we still 
can compute the K~ scattering length and which can 
be comparable with the above experimental value. 

Since this is a first lattice QCD evaluation of K^K~ 
scattering length, our primary goal is to show the ap- 
proach in a conceptually clean way. Checks for system- 
atic errors including the discretization error, finite vol- 
ume effect, and so on, are being investigated on larger 



MILC lattice ensembles: 24^^ x 64 and 28'^ 



96, 



48'^ X 144. Moreover, it will be important to repeat this 
calculation using more MILC lattice ensembles such that 
a continuum extrapolation can be performed. Further- 
more, the study of K'^ K~ scattering in the multi-channel 
scattering (e.g., tttt and KK two channels) [35| is un- 
der preparation, where we will explicitly consider the 
ao//o resonance [s^. We anticipate that these sophis- 
ticated lattice computations of strong K^K~ scattering 
length would be a specially important source for further 
progress toward a better understanding of some reactions 
like pp — > ppK^ K~ , and so on. 
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